Dynamical density matrix renormalization group study of photoexcited states in 

one-dimensional Mott insulators 
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One-dimensional Mott insulators exhibit giant nonlinear optical response. Based on the dynam- 
ical density matrix renormalization group method, photoexcited states and optical response in the 
insulators are studied as functions of the on-site and the nearest neighbor Coulomb interactions, U 
and V, respectively. We find that the lowest optically allowed and forbidden excited states across 
the Mott gap, which have odd and even parities, respectively, are degenerate for V/t 2 with t 
being the hopping integral of an electron between nearest neighbor sites. For V/t ^ 2, the bound 
states with odd and even parities occur and are not degenerate. The nature of the degeneracy and 
its effect on the optical response are examined. 

PACS numbers: 71.10.Fd, 71.35.-y, 72.80.Sk, 78.30.Am 
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In one-dimensional Mott insulators with strong on- 
site Coulomb interaction U, the ultrafast and gigantic 
nonlinear optical response have been observed [J 0- 
The third-order nonlinear susceptibilities in the typical 
one-dimensional Mott insulators of charge transfer type, 
such as cuprates and halogen-bridged Ni-halides, are a 
thousand times larger than those of band insulators or 
Peierls insulators [2] ■ In the previous theoretical investi- 
gations [1E|- it has been pointed out that the anomalous 
enhancement of the nonlinear susceptibilities in the one- 
dimensional Mott insulators originates in the large tran- 
sition dipole moment between an optically allowed state 
with odd parity and an optically forbidden (two-photon 
allowed) state with even parity due to the degeneracy of 
these states. 

In the large-t/ limit, only the charge degree of free- 
dom survives in the photoexcited states due to the spin- 
charge separation [j| , and the system is well described by 
an effective two particle model, called the holon-doublon 
model, which is composed of a holon representing one 
photoinduced empty site and a doublon representing one 
doubly occupied site with the attractive Coulomb inter- 
action — V between them. When V exceeds a critical 
value V = It, these paricles form an excitonic bound 
state 6]. Since these two particles cannot exchange each 
other, the photoexcited states with odd and even parities 
are degenerate regardless of the magnitude of V. By us- 
ing this model, the experimentally obtained optical spec- 
tra for the cuprates have been analyzed. However, the 
optical gap in the compounds is finite (~ 2 eV) and thus 
the holon-doublon model, which is obtained in the limit 
of large U, may be re-examined. Under a finite- U condi- 
tion, there is an intermediate process that these particles 
recombine, while this process is prohibited in the large-t/ 
limit. Since, in the finite-C/ condition, these particle can 
exchange their positions through the intermediate pro- 
cess, the eigenenergies of the odd- and even-parity states 



may be different. Therefore, it is indispensable to clarify 
the optical response in Mott insulators with finite U. 

The minimal model to describe photoexcitation in the 
Mott insulators is the single-band extended Hubbard 
model with U and V at half filling. By using the model 
with parameters for S^CuOs, the linear absorption and 
two photon absorption spectra have been calculated in 
the exact diagonalization method for a 12-site chain Q. 
To understand the degenerate feature in more detail, sys- 
tematic studies of the spectra are required for large size 
systems by changing the parameters widely. To achieve 
this purpose, we employ the dynamical density matrix 
renormalization group (DMRG) method @,H HE QH 
and calculate the dynamical correlation functions for 
odd- and even-parity states. The finite size scaling of 
the lowest-energy positions, which are derived from the 
corrections functions, is performed. We clarify a param- 
eter region where the degenerate photoexcitated states 
appear. 

The single-band extended Hubbard model in one- 
dimension is given by 
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where c\ a is the creation operator for an electron with 
spin a at site i, n% = + n^i, t is the hopping integral 
along the chain axis, U is the on-site Coulomb interac- 
tion, and V is the Coulomb interaction between nearest 
neighbor sites. We consider the half-filled case. We ex- 
amine the wide parameter region U/t = 5, 10, 20, and 
40 to see the nature of the excited states, although it is 
known that the realistic value of U/t for the cuprates is 
U/t ~ 10. 

We focus on the low energy region of the optical spectra 
with even and odd parities at zero temperature. The 
quantity Xa{&>) for an operator A is defined by 
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where L is the system size, |0 > is the ground state, Eq 
is its energy, and 7 is an infinitesimal positive number. 
Opeator A is either the paramagnetic current operator 
j = — it^2i a ( c l — H.c.) or the stress-tensor op- 
erator T = -tJ2i.a(4,a C i+l^ + H - C ') E2 • SinCe 3 and 

r have odd and even parities under the space inversion, 
respectively, Xjip) detects the odd-parity states while 
Xt(w) does the even-parity states. 

To calculate these quantities, we use the dynamical 
DMRG method 0, 0] , which is known to be powerful 
to study dynamical properties in one-dimensional corre- 
lated electron systems. The basic DMRG algorithm 0] 
for generating the ground state restricts the exponen- 
tially increasing Hilbert space to the relevant low energy 
DMRG basis using the ground state density matrix. To 
make informations about the photoexcited states reflect 
on the next DMRG step, the density matrix for a system 
block is expanded to a following form 0, E| : 



U/t=10 



U/t=20 



U/t=40 
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where ^2 a p a = 1 with a denoting each target state, the 
indices i and j run over all bases of a system block and an 
environment block, respectively. We adopt ip a =o — |0 >, 
ip a =i = >, and V Q=2 = (w + E -H + ij)- 1 Ai\0 >. 
The most relevant states to go onto the next DMRG step 
are constructed by the eigenstates of this enlarged ma- 
trix, and the spectral functions are directly given by the 
imaginary part of the inner product between ip\ and tp2 
after the sweep procedure is converged. We take po = 0.5, 
and pi = P2 = 0.25. It is noted that the precise vari- 
ational principle for the ground state wave function is 
necessarily violated whenever we treat the dynamical ef- 
fects by choosing po ^ 1. However, we can confirm that 
spectral shapes evaluated by DMRG are quantitatively 
consistent with results from the Lanczos diagonalization 
for small clusters. Since we empirically know that the 
numerical results are generally insensitive to how to de- 
termine the parameter set {p a }, we believe that conver- 
gency of the obtained ground state is considerably good. 

In the numerical simulation, an edge effect by the 
open boundary condition appears as a side peak located 
below a main excitonic peak. This indicates that a 
doubly occupied site tends to stay near the boundary. 
When we additionaly introduce an impurity potential as 
H = H + V{n\ + til), the side peak is removed. We 
have confirmed that the spectra obtained in the method 
is similar to those in the exact diagonalization method 
for small clusters with the peroidic boundary condition. 

M 

We first report the dynamical correlation functions, 
Xj(u) and Xt(u)i for several values of U/t and V/t in 
Fig. 1. The truncation number m is m = 180 in most 
cases of the simulation and sometimes m — 240 if nec- 
essary. The average discarded weight is about the order 
of 10~ 6 . In the largest-*//* case (U/t = 40), the spec- 
tra with odd and even parities are similar each other for 
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FIG. 1: C 1 -^) (solid lines) and x T (u) (broken lines) as func- 
tions of U and V. The system size L and the broadening 
factor 7 are taken to be L = 64 and "fL/t — 12.8. 



V/t < 3. With further increasing U/t at a fixed value of 
V/t, we can expect Xji^) — Xt(^) for any ui, since the 
model becomes equivarent to the holon-doublon model. 
In Figs. l(i)-l(l), a tiny bump structure is seen around 
lu = U — V, which is associated with a uniform exciton 
in the case of a dimerized spin background [l7^ . 

When the parameter U/t decreases at a fixed value 
of V/t, the odd-parity spectra, Xj( w )j does not change 
much. On the other hand, the spectral weight near the 
band edge in the even-parity spectra, % r (w), decreases, 
compared with that of Xj{w)- This is especially strong 

for V/t ^ 2 [see Fig. 1(d)] ^^|. As already mentioned in 
the introduction, the excitonic bound state is formed for 
V/t > 2 in the large-*/ limit. Even in finite-*/ cases, we 
expect that the bound state occurs as seen in the sharp 
peak structure of f° r V/t ^2 in Fig. 1. 

To evaluate the band edges of Xj( w ) an d Xt(w) more 
precisely, we carry out the scaling analysis in the case 
that U/t = 10. The band edge is defined by the lowest- 
energy peak position evaluated using a sufficiently small 
broadening factor 7. The scaling results are shown in 
Fig. 2. For 64-site systems, we take m = 300. For 
V/t < 2, it is found that the peak positions of Xj{<-°) 
and Xt(w) are almost degenerate at L — 64. The peak 
positions approach to the exact Bethe ansatz solution of 
the band edge at V/t = 0, which is shown by the bro- 
ken line in Fig. 2 |l9j. Since the extrapolated positions 
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FIG. 2: Finite size scaling for the lowest-energy peak positions 
of Xj{oj) and Xt(<^) for U/t — 10 . The filled circles denote 
the positions of Xi(u), and the cross symbols denote those 
of Xt(uj). The magnitude of error bars is smaller than that 
of the symbols. The band edge at V/t — obtained by the 
Bethe ansatz method is denoted by the horizontal broken line. 



in an d Xr(w) are insensitive to the magnitude of 

V/t for V/t < 2, we conclude the positions to be at the 
continuum band edges. On the other hand, the exci- 
tonic bound states are formed below the continum band 
edge with a certain gap for V/t > 2. In the case that 
V/t = 3, both parity states have the excitonic feature 
below the continuum band edge. However, the spectral 
weights around the bound states behave differently as 
seen in Fig. 1(d). Figure 3 shows the scaling result at 
U/t = 5. In this relatively small-[/ case, the critical 
value of V/t, where the degeneracy is lifted, still seems 
to exist near V/t = 2 but slightly smaller than that of 
the large-t/ limit. Consequently, we conclude that in the 
finite-?/ cases, the excitonic bound state formation for 
V/t ^ 2 lifts the degeneracy of the lowest-energy peaks 
in Xj( w ) an d Xr(w) seen for V/t < 2. 

The reason why the excitonic bound states are not 
degenerate in the case of finite U is interpreted as fol- 
lows: In the holon-doublon model, there is the local con- 
straint that double occupancy of a holon and a doublon 
is prohibited on the same site. On the other hand, in 
the case of finite U, there is an intermediate process that 
these particles recombine together. These particles can 
exchange their positions through the intermediate pro- 
cess. In such a case, nonzero matrix elements of the order 
of t 2 /U appear between photoexcited states before and 
after the exchange, leading to the lift of the degeneracy. 
For V/t £ 2, a probability of the exchange is large, since 
these particles approach each other to form the bound 
states. Therefore, the lift of the degeneracy becomes re- 
markable for V/t ^ 2. On the other hand, a probability of 
finding these particles nearby is expected to be small for 
V/t ^ 2, since the excitonic bound states are not formed. 
This indicates that the effect of the exchange is negli- 
gible. This idea is supported by the DMRG results for 
V/t<2, where the even- and odd-parity states are de- 
generate. We note that the behavior of the exciton for 



FIG. 3: Same as Fig. 2 but U/t = 5. 

V/t ^ 2 is similar to that of a Mott-Wannier type exci- 
ton |2(J , where an electron and a hole can exchange their 
positions freely. 

Returning to the extended Hubbard model, we clar- 
ify the mechanism discussed in the previous paragraph. 
To simplify the discussion, a strong-coupling perturba- 
tion theory from the large-t/ limit is applied to introduce 
an effective Hamiltonian. A unitary transformation to 
remove the electron motion from a singly (doubly) occu- 
pied state to a doubly (singly) occupied state leads to 

H = (U + J)^n ijT n a + (V - -J)^2mni + i 

i i 

i 

i 

where the three site hopping terms within the same Hub- 
bard bands are neglected, J = At 2 /{U — V), is 
the annihilation operator of the singly occupied state 
r\ l%a is the annihilation operator 

of the doubly occupied state r)i i(T = Ci t<T 7ii — a and Si is 
the spin operator. If we assume that all sites are occu- 
pied by one electron in the ground state, the photoex- 
cited states created by the current or the stress-tensor 



are written as (</>i ± 02 ) |0 > with 
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(Ei l(7 Vi+i,*Zi,<T, Ei, CT vlati+i,*)- We evaluate a propaga- 
tor matrix D(ui+ij) = t 2 < 0\<f>(oj+E -H+ij)- 1 $\0 >. 
Considering a sign factor by fermion exchange, we can 
regard — Im(Dn - D\2 — D 2 \ + D 22 )/L as XiW) an d 
-Im(£>n + D u + D 2 i + D 22 )/L as x t (cj). By use of 
the standard equation of motion method for the inverse 
Fourier transform of D, the optical spectra are finally 
given by 



XaM = -Im- 



4t 2 



ui + i-f — u!a - X A (uj + i-f)' 
where A is either j or r, u>j = ujq — J/2, lj t = loq 



(5) 
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J/2, w = U - V + 2 J, and < 0$ • S i+1 \0 > 
/ L = —3/4. The dynamical corrections are rewritten by 
Xj = (Y n - Y X2 )/2L, and X T = (Y u + Yu)/2L, where 
Y =< 0\8<f>(u + E Q - H + i7)~ 1 (5^ t |0 >j with I denot- 
ing the irreducible part of <f>. In a case that up and down 
spins on sites are alternately arranged by the interaction, 
J, 6ft is defined by 5ft = [2E il ,yS-i,.( ? )«»"li+i,»' ' 

For particular sites, i — i, i + 1, a transition process by 

is given by >= I • 0, T, IT, • • • >, 

which means that a holon-doublon pair extends spatially. 
Then, Xa corresponds to continuum band. Obviously, 
Xj and X T are degenerate at band edges according to 
their definitions. When the edges are located below uij 
and w r , the band edges of Xj and Xt are also degener- 
ate, since hybridization between the edge of Xa and a 
pole uja is negligible. We note that negligibly small hy- 
bridization corresponds to negligibly small probability of 
finding a holon and a doublon near by for V/t < 2 dis- 
cussed in the previous paragraph. When the excitonic 
bound states are separated from the edges, we can ap- 
ply the single-pole approximation to the calculation of 
X A , and obtain Xj(u) = X t (lu) = U 2 /[lu - (w + V)}, 
where we assume J2i < 0|jSf • 5^+2 10 > /L = 1/4. In this 
case, the energy difference between odd- and even-parity 
states is given by Aw = [u T + X T {tu T )\ — [uj + Xj(uij)] — 
J[V 2 - (J/2) 2 - U 2 ]/[V 2 - (J/2) 2 }. We note that the 
energy difference arises from the last term of the trans- 
formed Hamiltonian that describes exchange between 
a holon and a doublon. The values of Aw are evaluated 
to be Aw = 0.315 and 0.171 for (U/t, V/t) = (10,3) and 
(5,2.25), respectively. These values are comparable with 
the DMRG results for 64-site systems shown in Figs. 2 
and 3, which give Aw = 0.285 and 0.145, respectively. 
Since the exciton with even parity is situated near the 
continuum band due to Aw > 0, it mixes strongly the 
comtinuum band. Then, the excitonic spectral weight 
for the even-parity state is reduced. 

Among the related compounds, we focus on Sr 2 Cu03 
to make a comparision between the present results and 
experiments. Mapping the Zhang-Rice band onto the 
lower Hubbard band, the parameters in eq. (1) for 
Sr 2 Cu03 have been estimated to be U/t = 10 and 
V/t = 1.5 ~ 2.0 0,0]. The parameters correspond to 
those used in Fig. 1(c). We find that Xj( w ) representing 
the distribution of the odd-parity states reproduces the 



measured linear absorption spectrum yj, |2( . In the case 
that U/t = 10 and V/t = 2, the absorption band edges 
in Xj( w ) and Xr(^>) are found to be nearly degenerate. 
The distribution of the even-parity states for Sr2Cu03 
has been measured by the two-photon absorption (TPA) 
spectroscopy, and a finite energy difference (~ 0.3 eV) be- 
tween the peaks in the linear absorption and TPA spectra 
has been reported Q. The behavior of the band edges 
in the odd- and even-parity spectra is different from that 
of the peaks in the spectra. The calculation of the TPA 
spectrum in the present model will make it clear whether 
the peak position in the TPA spectrum is higher in energy 
than the band edge position in \ t (uj). Such a calculation 
is in progress. 

Among the halogen-bridged Ni-halides, the Br-based 
compounds show anomalous behaviors: the linear absorp- 
tion peak is very sharp but the excitonic effect is found to 
be small from photoconductivity measurements although 
the effect can not be neglected [22| . In the present model, 
the optical responses observed in the compounds have not 
been obtained. Since the optical gap A in the compounds 
has been reported to be comparable with the p — d trans- 
fer energy t p d |23|. it is not clear that the single-band 
model describes the electronic states in the compounds. 
It might be necessary to study the underlying electronic 
states more detail. 

In summary, we have presented the optical spectra 
with odd and even parities by changing widely a set of 
parameters in the one-dimensional single-band model by 
using the dynamical DMRG method. We have obtained 
three results: (1) The excitonic bound states are formed 
for V/t > 2 in finite-C/ cases. (2) When the excitons are 
unbound for V/t ^2, the photoexcited states with odd 
and even parities are degenerate. (3) When the exci- 
tonic bound states are formed, the eigenenergies of the 
odd and even states split because the empty site and 
the doubly occupied site can exchange each other in the 
finite- U cases. 
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